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1. Background 



This is the final report of the project "Study of 
Mark's Electrothermodynamic (ETD) Power Generator", spon- 
sored by the Department of Energy (DOE) . The primary ob- 
jective of this effort has been to provide a detailed 
theoretical analysis of a type of electrohydrodynamic (EHD) 
power generator proposed by the Marks Polarized Corporation. 
Once a sound analytical basis was developed, it was possible 
to ascertain the limitations of the proposed system and to 
seek improvements. 

The results of our first study are reported in detail 
in Ref. (1) and summarized more concisely in Ref. (2). 

That work indicated the originally proposed design would 
yield unacceptably poor performance and suggested a certain 
design revision. Further work was done to analyze such a 
design revision and is reported in the Appendix to Ref. (1). 
An improved simplified and shortened version of that work is 
reported in Ref. (3). Unfortunately, the new results indi- 
cated that the revised design, while showing somewhat im- 
proved performance, still would not appear sufficiently 
competitive in overall efficiency. 

In the first study, some minor effects of compressibili- 
ty of the medium were neglected. This was done in order to 
simplify the calculations. Also, it has been conjectured by 
Marks and others that it might be possible to improve greatly 
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the performance of EHD-generators which utilizes an ejector 
by employing a two-fluid cycle. Finally, performance gains 
that a diffuser might bring about appeared worthy of study. 

To this end, a second study for a two-fluid cycle was under- 
taken in which a diffuser is introduced just prior to the 
EHD generator and where full account is taken of the com- 
pressibility of both fluids at all locations. 

The principal earlier effort along the lines of our 
second study is Ref. (4), by Huberman, et. al , but our 
present work has gone beyond anything undertaken in Ref. 

(4). Although a two-fluid generator is analyzed in Ref. 

(4) , including the compressibility of both fluids, Huberman 
et. al, make no attempt to study the complete thermodynamic 
cycle. There are, moreover, a number of questionable 
assumptions which cast serious doubt in the validity of the 
results given in Ref. (4). 

The computational scheme that we have developed was 
reported in Ref. (5) . There a methodical derivation of 
the working formulas is given together with some preliminary 
results. Some further results are given in Ref. (6). 

The present report updates the material given in Refs 

(5) and (6) . In particular, we have found that instead of 
writing an energy balance across the condenser, it is more 
convenient to use a momentum balance as shown in Section II 

of this report. We have reproduced herein the complete analy- 
sis in its latest version. 
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A sample calculation for the two-fluid system is also 
included together v/ith the computer program listing. 

Figure 1.1 is a schematic of the generator configura- 
tion studied. In Marks' design, the centerline of the 
channel, which is shown in Fig. 1.1 as a straight line, is 
curved sc as to form a closed loop. 
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Fig.I.I Schematic of Generator Configuration 



Introduction 



2 . 



For purposes of analysis, the overall thermodynamic pro- 
cess in the two-fluid SHD generator may be represented by the 
schematic flow diagram shown in Fig. 2.1. 

Primary fluid enters the ejector at station 1, secondary 
fluid enters at station 2 and the mixture leaves at station 3. 

The mixture then passes through a diffuser, process 3-*-4 , 
the basic purpose of which is to reduce the kinetic energy of 
the flow and thereby decrease the associated friction losses 
downstream of the diffuser. On the other hand the diffuser 
introduces certain losses of its own so that some care is 
necessary to optimize the design in this respect. 

The gas mixture then passes through the electrical power 
section, process 4-*-5 , from which the gross electrical power is 
extracted. 

Finally, the gas mixture enters the condenser/separator 
at station 5. The primary fluid is separated from the gas 
mixture by condensation and leaves at station 6. The fluid at 
this point is assumed to be compressed liquid at a known stati 
pressure Pg and a known static temperature Tg . Kinetic 
energy at station 6 is regarded as negligible. The secondary 
fluid leaves the condenser/separator and enters the ejector at 
station 2. We assume, moreover, that 




( 2 . 1 ) 



and 




( 2 . 2 ) 



5 



Secondary Fluid 



3 

o 



o 

<D 

X 





6 



Fig. 2.1 SCHEMATIC FLOW DIAGRAM OF TWO - FLUID EHD CYCLE. 



The primary fluid which leaves the condenser/separator 
at station 6 is circulated back through the pump and through 
the boiler/superheater , process 7-*-l. 

It should also be explained that the primary fluid at 
station 1 is a condensible vapor of high molecular weight 
which contains fine liquid droplets. Unless specifically 
stated otherwise the quality z at this location is taken 
as 0.95. The important thermodynamic static properties at 
station 1 can be expressed in the following form. 

(2.3) 



h^ = h f + z h^ Enthalpy 

s^ = s^ - 1 -- z s^ Entropy 



(2.4) 



Here subscript f denotes the saturated liquid and subscript 
fg denotes the change between saturated liquid and saturated 
vapor. All properties are evaluated at a specified pressure 
; the corresponding saturation temperature T^ is then 
also known. 

The corresponding stagnation properties T g ^ and h g ^ 
are fixed when the stagnation pressure is specified. 

Of course the entropy has the same value s 1 at the stagna- 
tion condition as at the corresponding static state. Thus 

specification of P . and s. suffices to fix all other 

si 1 

stagnation properties. 

While it is convenient to start the analyses by stipula- 
ting the quality at station 1, it is also advisable later in 
the calculation to determine the corresponding quality at 
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station 4, the entrance to the electrical working section. 

It is essential that the fluid entering the working section 
contain finely dispersed liquid droplets of the proper size. 
These droplets carry the electrical charges that are the 
central feature of EHD power generation. We assume tentatively 
that acceptable droplet characteristics are obtained by 
maintaining the quality at the inlet to the working section 
in the range 0.92 to 0.98. These limits are only estimates. 

The thermodynamic analysis would be greatly simplified 
if both fluids were perfect gases with constant specific 
heats. In fact only the secondary fluid satisfies this re- 
quirement. Nevertheless, we can define an "equivalent perfect 
gas" which adequately approximates the essential thermodynamic 
properties of the primary vapor/liquid mixture. When carried 
out judiciously, this procedure simplifies the analysis at the 
cost of only a small loss in numerical accuracy 7 . 

The ratio of specific heats y of the "equivalent perfect 
gas" may be deduced from the important relation 



( p sl / p i> = < T sl /T l> 



Yj/fYx - 1) 



Solving this for 



y^ gives 



Y- 



1 



Z n(T sl /T 1 ) 

zHTp^Tp^T 



-1 



(2.5) 



( 2 . 6 ) 



Another fundamental property of the "equivalent perfect 
gas" is its specific heat C ^ which may be evaluated from 
the expression 
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(2.7) 



r = sl “r 
P 1 " (T sl - V 

where all quantities on the right are now known. 

It then follows from standard perfect gas relations that 
the gas constant and the molecular weight W, of the 

"equivalent perfect gas" are, respectively, 



R 



1 




( 2 . 8 ) 



and 




(2.9) 



where R = universal gas constant 
= 8317 joules/kg-mole °K 

Eqs. (2.6) through (2.9) complete the definition of the 
"equivalent perfect gas" which adequately simulates the essen- 
tial thermodynamic properties of the primary vapor/liquid 
mixture . 

It should also be explained that the secondary fluid at 
station 2 is predominantly a noncondensing gas of low molecu- 
lar weight primary vapor. 

It can be shown that the mass ratio co of condensible 
primary vapor to dry secondary gas at station 2 is given by 
the expression 



0) 




w 

l p 2 - W 1 



( 2 . 10 ) 
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where 



W v = true molecular weight of condensible vapor (not to 
be confused with the "equivalent molecular weight 
W^" considered earlier) 

W G = molecular weight of noncondensing gas 
P v (Tg) = vapor pressure of condensible primary fluid at known 
temperature Tg 

~ known static pressure at station 2 
The present analysis is restricted to conditions under 
which co is very small compared with unity and may be neglec- 
ted. Thus the fluid at station 2 may be treated as dry gas. 
This greatly simplifies the analysis. Eq. (2.10) is useful for 
verifying that this assumption is indeed satisfied in any par- 
ticular instance. 

It should be added that if it later becomes of interest 
to make calculations for circumstances under which parameter 
co is not negligible, the present analysis can in fact be 
generalized to include this effect. For the present, however, 
we prefer to deal only with the simpler situation in which co 
may be neglected. 
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3. Ejector 

The ejector receives a primary stream of high molecular 
weight gas at station 1, a secondary stream of low molecular 
weight gas at station 2 and discharges the resulting mixture 
at station 3. Static conditions at these three stations are 
designated by subscripts 1 , 2 , 3 . Stagnation conditions 

at the corresponding stations are designated by subscripts 
si , s2 and s3 . 

The following quantities are arbitrarily specified or 
known at stations 1 and/or 2, namely. 

Molecular weights: and W., 

Ratios of specific heats: y^ and y 2 

Stagnation pressure of primary jet: P gl 

Stagnation temperature of primary jet: T g ^ 

Static pressures: P^ = P, (= Pg) 

Static temperature of primary jet: T^ (= saturation tempera- 

ture of primary fluid at pressure P^) 

Static temperature of secondary jet: T 2 (= Tg) 

Mass flow ratio: m^/m^ = m^/(m^ + m 2 ) = x 

Velocity ratio: V^/V 2 = y 

We treat the fluids at stations 1, 2, 3, 4 and 5 as per- 
fect gases with constant specific heats. 

The velocity ratio of the ejector may be developed as 
follows: 



V 



V. 



1 



2 




M 



2 



1 



( 3.1) 
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We can now use this result and other standard perfect 
gas relations to complete the calculation of key thermody- 
namic properties at stations 1 and 2. It is convenient to 
arrange these calculations sequentially as follows : 



M i - Iff - 1 



= 






M, 



s 2 



= 1 + 



Y 1 R 1 T sl 

Y 2 R 2 t 2 T sl /T l 

Y 2 ‘ 1 „ 2 

^ 



s 2 



s 2 l 
5 2 J 

P 



t 2 (t s 2 /t 2 ) 

^2 

Jr, 

s 2 



T _.^2 “ X ) 



s 2 



= P 



s2\ 



2 P. 



(3.2) 

(3.3) 

(3.4) 
( 3.5) 

( 3.6) 
( 3.7) 



Before we can determine the actual conditions at station 
3, the outlet of the ejector, it is first necessary to analyze 
two other hypothetical cases as indicated schematically in 
Fig. 3 . 1 . In all three of these cases the flow is taken as 
steady and adiabatic. 

The first case, that shown in Fig. 3.1(a) , involves a 
hypothetical device which receives two gas streams at stations 
1 and 2 and discharges the resulting mixture at station x . In- 
let conditions at 1 and 2 are identical to those of the actual 
ejector. 
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Fig. 3.1 COMPARISON OF ACTUAL EJECTOR WITH IDEAL 
EJECTOR AND WITH HYPOTHETICAL MINIMUM 
LOSS DEVICE. 
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Recall that the adiabatic mixing of two different gases is 
an inherently irreversible process that always involves a cor- 
responding entropy increase. We stipulate that the first case 
is such that the only irreversibility which occurs is that asso- 
ciated with this mixing. We term this mass mixing to distin- 
guish it from another type of mixing considered below which we 
term momentum mixing. 

The second case, that shown in Fig. 3.1(b) , involves an 
ideal ejector which receives two gas streams at stations 1 and 
2 and discharges the resulting mixture at station y. Again, 
the inlet conditions at 1 and 2 are identical to those of the 
actual ejector. This case satisfies, among other relations, the 
idealized one-dimensional momentum equation for frictionless, 
constant area flow. It can be shown that in this case there is 
an overall entropy increase which includes not only the previ- 
ously mentioned effect of mass mixing but also a further increase 
associated with momentum mixing. 

The third case, that shown in Fig. 3..1(c) , represents the 
actual ejector itself which receives input streams at stations 
1 and 2 and which discharges the resulting mixture at station 3. 

All three of the above cases refer to constant area pro- 
cesses in the sense that 

A x = A y = A 3 = (A x + A 2 ) ( 3.8) 

Once the hypothetical operating conditions at stations x 
and y have been found, it then becomes possible to express the 
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corresponding actual conditions at station 3 . This is accom- 
plished through the use of an ejector effectiveness n E - 
This parameter is defined more explicitly in the later analysis. 
Values of n E must be estimated by reference to test data on 
ejectors (i.e.. Ref. 8, 9, 10) . 
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4. Adiabatic Mixing of Two Streams 

The process in the ejector is treated as the adiabatic 
mixing of two streams. The following relations can be shown 
to apply. 



R = 8315 Joule/kg'-rcole K 



R- = §- 
1 W n 



R, = 5_ 

2 W- 



R„ = R = R, = X R. + (1 - X) R„ 

X y 3 1 ^ 



Y 1 R 1 



pi " ( Yl - 1) 



Y 2 R 2 



P2 " 


(y 2 


- 1) 




C 


C 


■ c , 


= x C . + 


px 


py 


P 3 


Pi 
c - 


Y x = 


II 

>i 


u> 

11 


p 1 


<C p3 - R 3 


T 


T 


= T , 


= pi- [x 


sx 


sy 


s3 


C p3 


A x = 


II 

>1 

< 


A 3 = 


< A 1 + V 



P2 



'pi si 



- n 1 n 

p2 s2 



( 4-1) 
(4-2) 

( 4.3) 

(4-4) 

( 4.5) 

(4.6) 

(4. 7) 

(4.8) 

(4.9) 
(4.10) 



Notice that the quantities R , , y , T^ and A . all 

have the same values at stations x, y and 3. 

This fact may be used to simplify many of the subsequent 
expressions. It is always permissible to substitute subscript 
3 for subscripts x or y on any of these quantities. 
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5. Mass Flow 

It is convenient temporarily to omit station subscripts 
1, 2, 3, x, y and thereby develop certain needed mass flow 
relations in generalized form. Thus the mass flow across an 
arbitrary station may be written and developed as follows. 



where the auxiliary function f(M) is defined as follows. 



Eq. (5.2) can be applied specifically to stations 1, 2, 
x, y, 3 as follows 





P .A T /T 
s * s 



(5.1) 



This may be rewritten as 



P A 



it = — ^ f(M) 



(5.2) 





(Y + 1) 



f (M) 



(5.3) 




( 5. 4) 



m 



2 




(5.5) 
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(5.. 6) 



m = 
x 



P SX A 3 
^ R 3 T s 3 /y 3 



f (M ) 
x 



m = ^ ^ — f (M ) 

y , 3 — fP — itt y 



^3 T s3^3 



(5.7) 



m 3 = 



P s3 A 3 
l/ R 3 T s 3 /y 3 



f (m 3 ) 



(5.8) 



also 



ih. 



m. 



= x 



(5.9) 



m = ft = rti_ = iti, + rti 
x y 3 1 2 



(5.10) 



From Eqs . (5.4), (5.5), (5.9) and (5.10) we readily find 



that 



' P s2\I R 1 T sl h f(M 2> 



\ p s i/j R 2 T S 2 Y i f ^ M i^ 



(5.11) 



Then from Eq. (5.10) we infer also that 



(Ai/A 2 ) 

(1 + A 1 /A 2 ) 



(5.12) 



Since all quantities on the right side of Eq . (5.11) are 

now known, Eqs. (5.11) and (5.12) fix the area ratios A^/A 2 
and A^/A^ . 

From Eqs. (5.5), (5.6), (5.9), (5.10) and (3.8) we deduce 
further that 
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f <m 2 ) 



( 5 . 13 ) 



(1 - A /A-) 

f( V * a -V 




R 3 T s3 Y 2 
R 2 T s2 Y 3 



All quantities on the right side of Eq . (5.13) except 

P are now known. The method of determining P is ex- 
s x sx 

plained in a later section; refer to Eq . (5.15). Once P 

SX 

has been found, Eq . (5.13) fixes f(M x ) whereupon Eq . (5.3) 

fixes itself. The solution of Eq. (5.3) for M when 

f(M) is known involves an iterative procedure which is ex- 
plained later in this section. 

From Eqs . (5.6), (5.7), (5.8) and (5.10) we also find that 



P f(M ) = P f (M ) = P , f(M,) 
sx x sy y s3 3 



(5-14) 



Assuming P and M known, this relation along with Eq. ( 5 . 3 ) 

S X X 

fixes P when M is specified; it also fixes f(M_) and 
sy y 3 

when P s 3 is specified. The methods of calculating 
and P ^ are explained later; refer to Eqs. (7.8) and (8.14). 

Fig. 5.1 is a rough sketch of the function f(M) as defined 
by Eq. (5.3). This function vanishes at M = 0 and at M = « . 
The maximum value occurs at M = 1 and equals 



max 



y + 1 



(Y + 1 ) 
2(y - 1; 



(5 .15) 



It is clear from the figure that in the range 0<f< f , 
Eq. (5.3) has two positive real roots, one subsonic and the 
other supersonic. Usually the supersonic root must be discarded 
because it is not consistent with a further constraint imposed 
by the second law of thermodynamics. 
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Fig. 5.1 APPROXIMATE SKETCH OF FUNCTION f(M) 
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In order to solve Eq. (5.3) for M when values of f 

and y are specified, we utilize an iterative procedure based 

on Newton's method. Let M denote the nth trial value of 

n 

M . We then define a function F(M ) and its derivative 

n 

I 

F (M ) as follows, 
n 



F' (M ) 
n 



F(M n ) 



- 



= f 



1 + Y - ~ 1 M 2 
2 n 



(Y + 1) 

2"(y ~ 1) 



- M 



n 



(5.16) 



- 














(3 - 


(H-i) «„ 


1 + 


h - A* 2 1 


2 ( Y - 




V 2 / n J 





Y) 

T7 



- 1 (5.17) 



The (n + l)^t trial value of M may now be taken as 

(5.18) 



F(M ) 

M, . , . = M - , 

(n + 1) n 



F (M n ) 



The cycle of calculations defined by Eqs. (5.16) , (5.17) 
and (5-18) is repeated until the result for converges to 
a stable value at the desired level of accuracy. 

Thus 



M = Jiim M 

n 

n -*■ 00 



(5.19) 



It is readily apparent that the value of M so obtained 
satisfies Eq. (5.3) for the prescribed values of f and y • 
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6 . Entropy 

Consider the case shown in Fig. 3.1(a). The specific 
entropy s of the discharged gas may be expressed in the 

X 

form 




where and P x2 are the partial pressures of the two com- 

ponents of the gas mixture. Symbols P Q and T o denote the 
pressure and temperature of the ambient atmosphere. The en- 
tropy of each component is assigned the value zero at this 

reference state P , T 

o o 

Incidentally, the primary fluid, being condensible, might 
well exist only in the liquid state at the reference condition 
P Q , T q . Nevertheless, so long as it is in the gaseous state 
at station x , Eq. (6.1) may still be used. 

The partial pressures of the two components are propor- 
tional to the respective mol fractions. Hence 




X R 1 / P x\ 

X + (1 - X)R 2 [p / 

(! - x) r 2 /P x \ 
x R x + (1 - x)R 2 [P o J 




(1 - x)R 2 



R 



3 




(6 . 2 ) 



(6 .3) 
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Upon substituting Eqs. (6.2) and (6.3) into (6.1) and 



simplifying we may obtain the result in the form 




where 



(6.4) 



Sjjj - £n R^ - x £n (x R-^) - ( 1 - x) R 2 £n [(1 - x)R 2 ] 

(6.5) 



It can be shown that this quantity s^ represents the 
entropy increase associated with the mass mixing of the two 
different gases. 

Moreover, since the entropy of the stagnation state s 

X 

is by definition identical to that of the corresponding static 
state x , we may replace T and P in Eq. (6.4) by T 

XX S X 

and P , respectively. Thus we obtain 

b 2v 



X 



= s 



m 



+ C 



p3 





( 6 . 6 ) 



In previous studies of this series, we have considered 
only the special case in which the same fluid is employed for 
both the primary and secondary streams. In that case we have 



R 1 R 2 ^3 R 



c,=c 0 = c, = c 

pi p2 p3 p 



Y, = Yo = Yt = Y 



(6.7) 
( 6 - 8 ) 

( 6.9) 



Moreover, in this special case there is no mass mixing, so 
that we must set 
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0 



( 6 - 10 ) 



' s = 
m 



Consequently, Eq. (6.6) now simplifies to 



s x =0 + Cp tn I ^ I - R tn 



(6-U) 



Next reverting to the case of two different gases, let 
At denote a time interval during which unit mass crosses 
station x . During this same interval x units of mass cross 
station 1 and (1 - x) units of mass cross station 2. The 
entropy which leaves across station x must equal the entropy 
which enters across stations 1 and 2 plus the entropy increase 
caused by the mass mixing. Thus 



s = s m + x s, + (1 - x)s- 
x m 1 2 



( 6 . 12 ) 



where 



■ Si tn (W - R i fc 



(6.13) 



s 2 = C p2 * n \ T 



s2 



R 2 in \P 



s2 



(6.14) 



Upon eliminating s between Eqs . ( 6.6) and ( 6.12) , we 

find that s also cancels from the result. It is then a 
m 

simple matter to solve for the exit stagnation pressure in the 
form 



In 



>) - k ka - fe 



- X S 1 - (1 - X)s 2 



( 6.15) 
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Thus the sequence defined by Eqs. (6 .13), (6 >14) and 

(6 .15) new fixes P gx . Next referring back to Eq. (5 .13) , 
we can calculate the value of f (M ) ; finally, from Eqs. 

X 

(5 .16) , (5 .17) , ( S18) we can calculate M itself. Once 

X 

P , T = T _ and M are known , it is a simple matter to 
sx SX SJ X 

calculate the corresponding conditions P , T from the 

x x 

usual perfect gas relations. 
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7 . Momentum 

Consider the idealized one dimensional constant area 
ejector shown in Fig. 3.1(b). The momentum equation for this 
system may be written 



(P. + p. V. )k, + (P_ + p. V, ) A, 



(P + p V 2 )A (7.1) 

y y y y 



The mass conservation relation for this system can be 
expressed in the form 



pi A i 'i 

X 



p 2 a 2 v 2 

(1 - x) 



= p A V 

y y y 



(7.2) 



Let us divide the three terms of Eq. (7.1) by the corres- 
ponding three terms of Eq. ( 7 . 2 ). Notice that the areas can- 
cel out of the result. In this way we obtain 



<p i + °i v i 2) . „ , 

X + (1 - x) 



( p 2 + P 2 V 2 2 ) 
P 2 v 2 



(P + P V 2 ) 

y y y 

P v 

y y 



(7.3) 



Temporarily dropping the station subscript, we next deve- 
lop the typical term of Eq. (7.3) in generalized format as 
follows . 



(P + pv 2 ) 

oV 




(1 + yM 2 ) 

M /T 7T 
s' 




(1 + yM 2 ) 





T7“T 

— ) M 



(7.4) 



Observe that the pressure P cancels from the result. 
It is convenient to define the auxiliary function 
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g(M) - «- A + (IL^Ji) M 2 ( 7 .5) 

(1 + yM ) * 

With this notation Eq. (7.3) may be rewritten in the 

form 



X 7 R 1 T S1 (1 - x) A 2 T S2 = 1 /^3 T s3 

gTi-Tp"/ g(M 2 ) / y 2 g(M ) / y 3 



( 7.6) 



Rearranging gives 



g(M ) = / 



,/R- 



s3 



Y- 



g(M x ) 



-l 



si 



Y- 



(1 - x) 
g(M-) 



s2 



Y- 



( 7. 7) 



Since all quantities on the right are known, Eq. (7 .7) 

fixes g(M^) . Then itself can be found by inverting 

Eq. (7.5). Fortunately, an explicit solution is possible in 

2 

this case as Eq. (7.5) can be reduced to a quadratic in M 
The result is 



M 2 = (1 - 2yg 2 ) ± A - 2(y + l)g 2 
1 - y ( 1 2 Y g 2 ) 



(7. 3) 



The general character of the function g(M) is sketched 

in Fig. 7.1. Notice that g(M) vanishes at M = 0 , reaches 

its peak value g at M = 1 , and decreases toward the 

msx 

limit g as M-*» . It can be shown that 

^ CO 



CT 

•"max 



1 

/2 ( y + 1) 



g 



CO 



l /(y -""T7 

Y / 2 



(7.9) 



(7.10) 
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Fig. 71 APPROXIMATE SKETCH OF FUNCTION g(M) 
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Fig. 7.1 discloses the existence of a supersonic root over 
the range g oo < g < g max . This root should not be accepted until 
a check calculation is made to determine whether the result is 
consistent with the second law of thermodynamics. The subsonic 
root is found by retaining only the negative sign before the 
radical in Eq. ( 7 . 8 ). 

Once M has been found from Eq. ( 7 . 8 ) , the corresponding 

3) and (5.14) . This 
y and sv . 



value of P may be found from Eqs. (5 
suffices to fix all properties at states 



29 



8. Availability 

The entropies s, and s 9 of the two input streams 
were defined earlier in Eqs. (6.13) and (6.14). The corres- 
ponding steady flow availability functions with respect to an 



ambient atmosphere 
be written 


at pressure 


P 

o 


and temperature 


T q may 


h 


= C , (T . - 
pi si 


V 


- T s . 
o 1 


(8.1) 




= C _ (T _ - 
p2 s2 


V 


- T s. 
o 1 


( 8.2) 



Consequently the total available energy entering the sys- 
tem becomes 



i|J m = x 'i'l + U - x )^ 2 



( 8.3) 



The respective entropies of the streams leaving at stations 
x , y and 3 are 



s = 
x 



s + C , 
m p3 




( 8.4) 



s 



y 



s + C 
m p3 



in 




R, Jin 




(8.5) 



s 



3 



s 

m 



+ C 



p3 





( 8 - 6 ) 



where s^ has been previously defined in Eq. (6.5). 

The corresponding availabilities of the streams leaving 
at stations x, y and 3 may therefore be written 



30 



: P3 


(T s3 " 


V 


- T 

o 


s 

X 


(8.7) 


: p3 


(T s3 - 


V 


- T 

o 


s 

y 


(8. 8) 


p3 


(T s3 - 


V 


- T 

o 


s 3 


(8.9) 



Now consider the losses of availability listed below. 
By substituting the above expressions for the quantities on 
the left and simplifying we readily obtain the expressions 
shown on the right. Thus 



(4* ~ 40 = T s 

m x o m 



( 8 . 10 ) 



( 4 > - 'l> ) 

X V 



R, T Jin . 

3 o \ P 



sx 



( 8 . 11 ) 



sy, 



<*X - 



R T In ( 

3 ° \ p s3, 



( 8 . 12 ) 



These results are very significant. Eq. (8.10) defines 
the loss of available energy caused by the mass mixing of the 
two different gases. Eq. (8.11) defines the further loss of 
available energy caused by the momentum mixing in an idealized 
ejector. Eq. (8.12) defines the augmented loss of available 
energy caused by the momentum mixing in the real ejector. 

It is useful to postulate an empirical relation between 
these last two losses which is shown below on the left. The 
expression on the right then follows from Eqs. (9.11) and 
(8 .12) . Thus 



<*x - V (p sy /P sx> 

<*x - *3> ' n E ‘ ^ ''W'W 



( ai3) 
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We term n E the ejector effectiveness. Its value, 
which is always less than unity, must be estimated from appro- 
priate test data on ejectors. 

If we treat n E as known, Eq. (8.13) fixes P ^ . The 
solution is simply 




We can next find f(M^) from Eq. (5.14) and from 

Eqs. (5.16), (5.17) and (5.18). This suffices to fix all pro- 
perties at states 3 and s3 thereby completing the analysis of 
the ejector. 
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Diffuser 



It is convenient to analyze the diffuser in terms of the 

concept of availability in steady flow. However, we wish to 

deal with the characteristics of the diffuser itself and these 

cannot depend on the arbitrary values P q , T q which happen 

to characterize the condition of the ambient atmosphere. We 

can achieve our objective by defining the availability with 

respect to a reference state P^ / T^ which characterizes 

the diffuser rather than with respect to the ambient state 

P , T . 
o ' o 

Under these circumstances we may write the appropriate 
availabilities at diffuser inlet and outlet as follows 



These equations make use of the fact that T g4 = T ^ • 
Also, states s3 and 3 are at the same entropy so that 






From Eqs. (9.1) and (9.3) we conclude that 
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This result shows that the available energy at the 
diffuser inlet is simply the inlet kinetic energy itself. 

The loss of availability through the diffuser may be 
found by subtracting Eq. (9.2) from (9.1) . The result is 



( ^s4 ' W = W = + R, T, Jin f ^ 



3 3 



s4 



We now define the diffuser effectiveness as 



s4 



'P 



= n, 



s 3- 



so that 



/ K 3 ~ ♦gA . 

\ / »S3 

This reduces readily to 



(1 - n D ) 







Solving for the pressure ratio gives 




e 



(l - n D ) 




This is the result required. It fixes P g4 when P 
and n D are specified. 



(9.5) 
( 9.6) 

( 9.7) 
( 9.8) 

( 9.9) 
s 3 ' 
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Assuming the area ratio A^/A^ to be arbitrarily spec- 
ified, we can find from the continuity relation. Refer- 

ring to Eq. (5.8) we may write 






P , A 0 
s3 3 



T s3 



7T- 



f(M 3 ) = 



P s4 A 4 
/ S 3 T s3 /Y 3 



f(M 4 ) = m 4 



(9.10) 



Consequently 



f(M 4 ) 

Eq. (9.11) fixes 
way from Eqs. (5.16) , 
to fix all properties 




f (M 4 ) . Then M 4 

(5.17) and (3.18) . 
at states 4 and 



(9.11) 

follows in the usual 
The result suffices 
s4 . Thus 





1 + 




( 9 . 12 ) 



= T ,/(T ,/T.) 
s3' s3 / 4 



(9.13) 



Y 3 




P 4 " 



(9.15) 



The configurations studied in earlier work did not incor- 
porate a diffuser. For the purpose of comparing the results 
of the present analysis with corresponding earlier work, it is 
desirable to be able to eliminate the effects of the diffuser 
in certain cases. This can be accomplished in the present 
analysis simply bv setting 
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1 



(9.16) 



n 



D 



and 




- 1 



(9.17) 



It is easy to confirm that under these circumstances 



s4 


= T s3 


(9-18) 


s4 


= P s3 


(9.19) 


M 4 


= M 3 


( 9. 20) 


A 4 


= etc. 


( 9.21) 



and all diffuser effects disappear. 
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10 . 



Electrical Power Section 



The power conversion process in an EHD duct is usually 
treated on the basis either of constant area or of constant 
static state. Inasmuch as electrical power output per unit 
mass is small, the numerical differences between the results 
computed by these two methods is negligible. For definite- 
ness in this analysis, however, we assume constant area. 

The electrical power that can be obtained from an EHD 
duct of constant area, negligible change of density and optimum 
length can be estimated from the one dimensional version of 
Poisson's equation which governs the electrical field. The 
solution is well known and will not be derived here; a detailed 
derivation may be found in Ref. (3). The essential result may 
be written in the form 



P 





V 



4 



where 

P g = gross electric power output, watts 

e = permittivity of medium 
— 12 

= 3.854 x 10 farad/m (for any gas) 

E, q = dielectric strength of medium, volts/m 
A 4 = area of duct (constant) 

= velocity 



( 10 . 1 ) 



According to the test data analyzed in Ref. 4, the dielec- 



tric strength is well approximated by the expression 




( 10 . 2 ) 



where C b3 is a characteristic constant of the medium. The 
data show that Eq. (10.2) applies to air or steam up to about 
10 atmospheres pressure with 



In this study we assume that Eq. (10.2) can be extrapolated 
up to about 100 atmospheres pressure. We also assume that an 
expression of the same form applies to other media besides air 
or steam but that each medium has its own characteristic value 
of the breakdown constant C Q . 

In the two-fluid system, the separate values of the 
breakdown constants and are usually known, but the 

breakdown constant °f the resulting mixture is seldom 

known. In the absence of adequate test data bearing on this 
point, we tentatively assume that can be estimated from 

the hypothetical relation 



Upon substituting Eq. (10.2) into (10.1)and dividing 
through by the mass flow rate, we find the electrical work out- 
put per unit mass of fluid in the form: 



C b3 = 9.49 x 10 3 m 2 u K/cmb 



(10.3) 
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By use of standard perfect gas relations, Eq. (10.5) ca n 
be changed to a more useful form. In this connection set 



(T 



s4 



T ,-) 
s5' 



AT 

s 



( 10 . 6 ) 



0 C P 

8 = — ^ — 2 . (io.7) 

T o 

2 

where P Q = ambient pressure, N/m 

T q = ambient temperature, °K 

Notice that 8 represents a dimensionless version of 
the electrical breakdown constant. 

Eq. 10.5 now yields the important result 



1 




The corresponding gross electrical work output per unit 
mass of fluid is then simply 

W * = C_ , AT (10.9) 

e P3 s 

Unfortunately, the known value of e and the typical 
experimentally measured values of are very small. Con- 

sequently, the dimensionless breakdown constant 8 is ex- 
tremely small compared to unity. Hence the drop in stagnation 
temperature AT e and the gross electrical work per unit mass 
W* e are also typically very small. This is a very basic and 
serious limitation on the performance that can be obtained from 
an EHD generator. 

Treating the process 4 -*-5 as isentropic, we may write 
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( 10 . 10 ) 



^3 ' 1} 



s5 



= P 



s4 \ T 



's5_ 

s4, 



By analogy with Eq. (5.8) we write the continuity rela- 
tion as 



P . A. P _ A c 

*4 = ^ 4 f(M 4 ) = 55 5 f(M 5 ) = * 5 



( 10 . 11 ) 



/*3 T s 4 /y 3 



v^3 T s5 /Y 3 



It then follows that 



^ ‘VS fe) ‘ ‘V 



( 10 . 12 ) 



This fixes f(M^) whereupon .M^ may be found in the 
usual way from Eqs. (5.16), (5.17) and (5.18). However, con- 
vergence should now be very rapid since we may take as a 
first approximation 



M_ ~ M 
5 4 



(10. 13) 



This solution now suffices to fix all properties at 
states 4 and s4 . 
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11. Condenser/Separator 

The primary fluid which leaves the ccndenser/separator 
at station 6 is assumed to be compressed liquid at known 
static pressure P g = P 2 = P x and at known static tempera- 
ture Tg = T 2 * Th e kinetic energy of the liquid at station 
6 is assumed to be negligible. 

Temperature T g must be equal to or greater than the 
ambient temperature T q in order to satisfy the require- 
ments of a heat balance on the condenser/separatcr . No 
attempt is made in this analysis to formulate this actual 
heat balance; instead T g is simply treated as a given or 
known quantity. However, it may be of interest later to 
investigate the effects of varying T g on the overall per- 
formance of the EHD system. 

The secondary fluid which leaves the condenser/separator at 
station z consists primarily of a low molecular weight nonconden- 
sing gas. It also contains some secondary vapor, but the mass 
fraction of this vapor is treated as negligible in this analysis. 
Thus the secondary fluid at station z is treated as dry gas. 

The thermodynamic process in the condenser/separator can be 
most conveniently analyzed and modelled by employing the momentum 
theorem. We treat the condenser/separator essentially as a one 
dimensional, constant area device in steady flew. By equating 
the net applied pressure and drag forces to the change of 
momentum flux in the direction of flew, we obtain the following 
simple relation, namely. 
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D = (1 - x)m,V 
3 z 



( 11 . 1 ) 



(p 5 - p z )a 5 




The effective drag force D in this equation represents two 
irreversible momentum loss effects. One of these is the loss of 
momentum produced when the condensible vapor in the incoming flow 
is condensed and brought to a halt on the stationary surfaces of 
the condenser/separator. The other is the additional friction 
drag force; the latter is conveniently expressed in terms of an 
empirical friction coefficient c f . Thus we may write 



Next we substitute Eq. (11.2) into 11.1) , bring all terms 
to one side of the equation, regroup terms, and divide through by 
the quantity p^A^ . We also make use of the following auxiliary 
relationships, namely. 




( 11 . 2 ) 




(11.3) 



and 




(11.4) 



In this way the following result is finally obtained - 



E (x) 




(1 - x) M 




= a (11.5) 
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The quantity E(x) serves here merely as a convenient abbre- 
viation for the sum of terms shown. It may be said to represent 
an excess of momentum available over that required to satisfy the 
various losses that occur across the condenser/separator. In any 
case, for assigned values of all other parameters, we seek by 
successive trials a value of x such that E(x) = 0 . Thus Eq. 
(11.5) fixes the mass fraction x in a manner which is consistent 
with the other specified input parameters. 

Of course, if the various fixed input parameters are not 
suitably chosen, Eq. (11.5) cannot be satisfied for any trial 
value of x in the range 0 < x < 1 , and no solution is physi- 
cally possible under these conditions. 
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12-. Cycle Efficiency 



Once a solution has been found such that E (x) = 0 , it 
becomes a s traightfoward matter to calculate the correspond- 
ing overall cycle efficiency. 

Because of the pressure drop through the boiler/super- 
heater , the pressure at boiler inlet is slightly 

higher than the stagnation pressure P g ^ of the primary 
fluid at the ejector inlet. We assume that the ratio 
is a specified constant. Hence with P g ^ specified, P^ 
is also known. Also the pressures Pg = P^ are known. 



Consider the ideal pump work w^ done per unit mass of 
primary fluid. For the present consider the hypothetical case 
of a reversible pump. Also note that density changes across 
the pump are negligible. Hence we may write 



* 



w 



P 




( 12 .. 1 ) 



The ideal gross electrical work output per unit mass of 

mixed fluid has earlier been established in Eq.(10.5). It is 

* 

denoted by symbol w^ . 

Let us now denote the net useful electrical work output 



per unit mass of primary fluid by symbol w net • It may be 
* * 

related to w g and w^ in the following way 



w 



net 




Joule/kg 



( 12 . 2 ) 



where n , the excitation efficiency, allows for the small 

X 

electrical power expended to excite the system and where n , 
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the pump efficiency, allows for irreversibilities in the pump. 
The heat input in the boiler/superheater per unit mass 
of primary fluid may be written 

★ 

>Jin. * (h sl - V ' W 2 (12 ' 3 > 

P 

where hg is the enthalpy of the primary liquid at conditions 
P g , Tg as listed in suitable tables of properties. 

The overall thermal efficiency of the cycle can now be 
calculated from the simple formula 



n 



c 



w 



net 

< 5 * 



m 



(12.4) 



Both n and w are useful parameters which charac- 

c ne t 

terize the overall thermodynamic performance of the EHD sys- 
tem. 

Eq. (12.2) shows clearly how the ejector serves as a 
kind of amplifier which increases the electrical work output 
per unit mass of primary fluid. This effect is shown by the 
presence of the parameter x in the denominator of the first 
term. Notice the beneficial effect of a lew value of x on 
parameters w ^ and n c • Unfortunately, the value of x 
cannot be stipulated independently in advance; it is fixed 
by the other specified input parameters as explained in the 
previous section. 

★ 

Notice that owing to the fact that w g is normally 
very small, it is possible in some circumstances for w net 
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to take on values which are actually negative. When this 
happens it means that electrical power output is insufficient 
to drive the pump. 

The relationships developed in this and the preceding 
sections make it possible to carry out systematic parametric 
studies of various one-fluid and two-fluid EHD cycles, with 
realistic allowances for the various losses that occur. Such 
studies can establish optimum design parameters and perfor- 
mance limits under various circumstances. These results in 
turn can finally permit informed conclusions to be drawn con- 
cerning the ultimate feasbility of this general type of EHD 
power generator. 
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13. Calculation Sequence 

In this section we summarize the various equations pre- 
viously derived in the approximate order in which they would 
be used in the calculation of system performance. Also listed 
are the initial input parameters whose values must be specified 
in order to start the calculation and various further input 
parameters whose values must be specified in order to proceed 
with various successive stages of the calculation. 

Initial Input Data : p i = p 2 = P 6 ' T 1 ' P sl ' T sl ' ^ h sl “ h l^ ' 



R = 8315 Joule/kg°K 




( 13 . 1 ) 




(13 .2) 




(13 .3) 




1 



(13 .4) 



Further Input Data : W , W G = W 2 , T 2 = Tg , P v (Tg) 




(13 .5) 



Verify that cocci 
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Verify that t^l 




Further Input Data: x (trial value) 




R 



xR^ + 

Y 2 R 2 
(y 2 - 



(1 



1) 



- x) R 2 



(1 - x) C . 

P2 



( 13 . 6 ) 

( 13 . 7 ) 

( 13 . 8 ) 

( 13 . 9 ) 

( 13 . 10 ) 

( 13 . 11 ) 
( 1 . 3 . 12 ) 

(1 3 . 13 ) 
( 1 . 3 . 14 ) 

( 13 . 15 ) 

( 13 . 16 ) 
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(13.17) 



Y3 = (C p3 * 3 r 3> 



T s 3 “ lx C pl T sl + (1 - x > C p2 T s2 J 



f (M x ) 



= M, 



Y 1 “ ^ -7 

1 + — 2- w r 



(Y x + 1) 
2(7! - 1) 



f(M 2 ) = M 2 



1 + 



y 2 - 1 



M. 



( Y 2 + 1) 
2(Y 2 - 1) 



A. 



x / P s2\ Pi T sl ^2 f(M 2> 

, 2 - (1 - x) \P gi ; v/R 2 T s2 Yl I(M x ) 



A, 



(a 1 /a 2 ) 



a 3 (l + a 1 /a 2 ) 



Further Input Data: P , T 

o o 



-x - Cpl t» (y - *1 (5f 



2 p2 



i /p . 

F f j - s 2 la \t 



s, = C_, In ( - 3 2 2n [ _ s2 



in (W’^! c P 3to (W' xsi ' a ' x) P 



in (P sx /P 0 ) 



P = P e 
sx o 



(13.13) 

(1 3.19) 

(1 3.20) 
(1 3.21) 
(1 3.22) 

(1 3.2 3) 
(1 3.24) 
(1 3.25) 

(13.26) 



49 



(13.27) 



(1 s2\ i|^3 T s3 Y 2 

f< V - £ a - xi fc) \] r 2 t s2 Y3 £ ' m 2 > 



(Y 3 + 1) 

2 \ 2 '^Y 



'max ly-, + 1 



Verify that f(M ) = £<£ 

x max 

For first approximation set M = 

X z 

Iterate using Eqs. (13.31), (13.32), (1,3. 33). 



F(M n ) = f 






(y 3 + 1) 
2(y 3 - l) 



- M 



n 



F' (M n ) = f 



fr: 



+ i) ( 

— y 






(3 - Y 3 ) 
2(Y 3 ~ T) 



F(M ) 

M (n + 1) M n " F' (M n ) 



Iteration converges to fix M, 



x 



M 1 \ Y 1 " 1 2 
g( M ) = L—. 1 + (-L- ) M, Z 

(1 + y M ± Z ) V 2 1 



g(M 2 ) = 



M- 



Y 2 ^ .. 2 



(l + y 2 m 2 ) 



2 ^ \| 1 + ^ 2 ^ M 2 



(13 .2 8) 

(13 .29) 
U 3.30) 



(13.31) 



(13.32) 



(13.33) 



(13.34) 



(13.35) 
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v-1 



g(w y ) 



= g = 



R 3 T s3 



x 



Y-, 1 g(M, ) / Y 



R-, T , , . > 

1 si ( 1 - x) 



g(M 9 ) 



R T 



s2 



(13.36) 



M 



2 = 



(1 - 2y 3 g ) - y 1 - 2 (y 3 + 1) g 

l - y, (l - 2 Y , g 2 ) 



(13.37) 



Eq. (13.37) applies provided that g ^ g^ = — \ 



(Y, ~ 1) 



(13-33) 



If g = g OT the solution reduces to 



( Y -? - 1) 



H y = 



2y- 



(13. 39) 



Once M y is known from Eq. (13.37) or (13.39) 



f(M ) 

y 



(y 3 “ J-) 2 

= M y l 1 + -h M y 



( Y 3 + 1) 
2 ( Y 3 " 1) 



(] 3.40) 



f( M ) 

> = p ±_ 

sy sx f(My) 



(13.41) 



Further Input Data: 



(P 



l/n. 



= ? 



s3 sx 



sy i 

SXJ 



f(M 3 ) - ? 



SX 



s3 , 



f(M x ) 



Verify that f(M 3 ) = f<f 



Max 



(13.42) 

(13.43) 
(13. 44) 
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(13.45) 



For a first approximation take = M 

Iterate using Eqs. (13.31); (13.32), ( 13.33). 

Iteration converges to fix . 




1 + 

. 




M. 



2 



T s3 /<T s3 /T 3> 






1 ) 



" P s3/ (P s3 /P 3> 



Further Input Data : n D 




P 



s4 





f (M 4 ) = 




Verify that f(M.) = f<f 

4 mcix 

For a first approximation take 

Iterate using Eqs. (I 3 . 3 I), (13.32) 

Iteration converges to fix M . 




(13.33) . 




(13.46) 

(13.47) 

(13 . 48) 

(13.49) 

(13.50) 

(13.51) 

(13.52) 

(13.53) 



(13.54) 



52 



(13.55) 



' s4 



. i ♦ fl24i) 



M y 



T 4 - T s4 /(T s4 /T 4 ) 



(13.56) 



P T 3 ” 1 ) 

F s4 X s4 J 



P 4 " P s4 /(P s4 /P 4 ) 



(1.3.57) 

(13.58) 
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Further Input Data : e = 8.854 x 10 farad/m, Cg^ , Cg 3 > 

3 9 

(C n = 9.5 x 10 m“ °K/cmb for air or steam) 

U 



C B3 R 3 R 1 C B1 + (1 “ x) R 2 C B2^ 



(1 3. 59) 



3 = 



£ C B3 P o 



(1 3. 60) 



AT 



' s4 



= 3 



3 

, AT 

_s5 = (1 - s- 1 ) 

■ „ s4 

s4 



^3 - l \(*s4\ro 
2 y 



2 / 



s 4 j 



Y3 ~ 1 9 

1 + (^-2— )M 4 



-i/(y 3 



- 1 ) 

(1 3. 61) 

(1 a 62) 



s5 



T „ 

= P . (=^2) 



y 3 /(y 3 - i) 



s4 T 



(13. 63) 



s4 
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f(M 5 ) = f = 



(13-64) 




Verify that f(M,-) = 

For a first approximation take 

Iterate using Eqs . (1 3-31), (1 3-32), (1 3.33) 

Iteration converges to fix . 




Further Input Data : c f 

p ^ ‘F 2 

E(X) = (1 “ -~) + (1 - X - — )Y 3 M, 

p 5 2 J ^ 



- (1 - x)y 3 M 5 M, 



y 2 r 2 T 






Y -3 *7 T 



(13-65) 

(13.66) 



(13-67) 

(13.68) 



(13-69) 

(13.70) 

(13 .71) 
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Iterate from Eq . (13.14) until x converges. 



Further Input Data: (P sl /P g7 ) / ° 6 * (h gl - h g ) , n x , n p 



P _ = P ,/(P ,/P -) 
s 7 si si s / 



w = C , At 
e p3 s 



(13.72) 

(13.73) 



w 



(P - P ) 
* = ' s7 V 



(13.74) 



w 



n w 
x e 



w 



ne t 



x n. 



RESULT (13.75) 



w 



[ in (h sl " n 6 } " tp 



(13.76) 



w 



n c q 



net 



m 



RESULT (13.77) 
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14. OPTIMUM VELOCITY RATIO 



The calculation sequence summarized in the preceding sec- 
tion reveals that a comple cycle calculation requires us to 
specify numerical values for about thirty input parameters. 

One of these is the velocity ratio (V^/^) = y . The calcula- 
tion procedure outlined above then permits us to find the cor- 
responding value of the mass flow ratio x , assuming that 
there exists a value of x which satisfies the governing 
thermodynamic relations for the specified values of the other 
input parameters. 

If we now change the value of y but hold all other in- 
put parameters fixed, we can repeat the above procedure and 
find a corresponding new value of x , if such a value exists. 
It is evident, therefore, that x becomes some definite 
function of y as long as we remain in the domain where a 
real solution exists. Of course all other dependent variables 
of the cycle are also functions of y . In particular, the 
overall cycle efficiency ri is some definite function of y 
Moreover, for prescribed values of the other input parameters, 
there will be some definite value of y , let us designate it 
as the optimum value Y 0 p t / which yields the greatest value 
of h c that is possible under the specified input conditions. 
Our basic purpose therefore is to determine for any prescribed 
values of the other input parameters, the value of Y 0 p t an< ^ 
corresponding value of n c ) max . Of course the- values of x opt 

and other dependent variables are also of interest. 



the 
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Since a working computer program is now available to es- 
tablish values of y ^ , n ) , x and so on for pre- 

^opt c max opt ^ 

scribed values of the other input parameters, we can proceed to 
study the effects of changes in these other parameters and to 
search systematically for such values of the most important 
input parameters as will yield the best overall performance 
of the system. 
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15. Conclusions and Recommendations 



Based on the work reported here and in Ref. (6) , the 
following conclusions are evident: 

1) Ejector performance is substantially improved 
through the use of a high molecular weight pri- 
mary fluid (i.e., above 100) and a low molecular 
weight secondary (i.e., below 30). The mixing 
penalty which ensues is offset by the substantial 
decrease in ejector losses. 

2) The breakdown strength limitation is severe even 
at higher pressures. Since the breakdown strength 
of a mixture can often be smaller than that of a 
single fluid, allowance for this must be made when 
considering fluids to be used in the ejector. In 
Ref. (6) it is stated that in mercury /hydrogen an 
improvement by a factor of 5 in the breakdown 
strength is required for a competitive thermodynamic 
efficiency. 

3) Consideration of the breakdown strength of mixtures 
together with other factors seem to indicate that 
the sequence of components shown in Fig. 1.1 might 
not be the best possible. In particular, serious 
thought should be given to an arrangement where the 
electrical section precedes the ejector. Here, 
mostly secondary fluid is present and the breakdown 
strength will be at its highest value. 
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The termination of the contract with DoE has left some 



questions abcut the two-fluid cycle unanswered. With the 
analysis presented in this report it should be possible to 
carry out further thermodynamic studies with various fluid 
combinations. Also, more fundamental knowledge on the 
breakdown strength of dielectric fluids will be necessary 
before the prospects for successful EHD power generation 
can be fully evaluated. 
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16. Sample Output 

A sample output is shown in the following pages. The 
case is for mercury /hydrogen expanding from 1500 psia to 
700 psia. The value of the breakdown constant, Cg, is 8000. 
The component loss parameters have numerical values corres- 
ponding to CASE B in P.ef. 6. 
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FLUID PROPERTIES 



1 


PS1= 1500.0 PSTA 


p 


TR1=s 1*45.0 UFn 


P 




3 


Pl= 700.0 PSIA 




4 


n- 1770.0 npn 


R 




cr 


< MS 1 -HI > = 13.50 


RTIt/LB 


A 


(H81 -HA )= 1 70 70 


RTIj/l ft 


7 


RH06= 0 . 7500E 03 


LR/FT3 


8 


r R 1- O.ROOOF 04 


M7 


OFF K/PMR 


9 


NFLIJID= 2 




to 


1*1(3 = O,-<0 TAP oi 






11 


GAMMA2= 1.4000 




12 


r: ft 7s o.rooof 04 


M2 


npn K/rMft 


SYSTEM PARAMETERS 


\ 










1 


ET AE= 0.85 2 


ET AD= 0 ♦ 85 


3 


A3/A4=0 . 250 4 


CF = 


c 

oc 

c 


5 


ETAX= 0.98 6 


ETAP- 0.8*5 


7 


PS1/F*S7= 1*00 






COMPUTATIONAL FORMAT 












1 


DELY= 0.10 2 


YMIN= t.00 


7 


YMAX = 5 . 00 







UV= 0.2058E 03 ( EFF MOL UT VAPOR) 0AMMA1= 1 .IO 1 A Y:.»r»= 0.10A4 
ENTER CONTROL CODE 



Y 


X 


A1/A3 


OOFPAI. 1 
FFF 


RF! Al TUP 
EFF 


M3 


M5 


5.00 


0.4375 


0.0015 


0.0 228 


0.9285 


0 . 0785 


O . 0071 


4.90 


0.4531 


0.0017 


0.0213 


0.2175 


0 . 0?95 


0 . 0074 


4.80 


0.4375 


0 . 0 0 1 6 


0. 0228 


0.2285 


O , 099 A 


0 . 0074 


4.70 


0.4375 


0.0016 


0.0228 


0. ?or5 


0.0303 


0 * 00 7 A 


4.60 


0.4277 


0.0016 


0.0238 


0.2384 


0.0707 


O > 0077 


4.50 


0.4063 


0.0015 


0. 0261 


0.2A20 


0.0308 


0.0077 


4.40 


0.4219 


0.0016 


0.0244 


0.7446 


0.0319 


0 . 0080 


4.30 


0.4063 


0.0016 


0.0261 


0.2620 


0 . 03°9 


0 . 0080 


4.20 


0.4063 


0.0016 


0 . 0261 


0 . 2 A90 


O , 0379 


O , OOR7 


4.10 


0.3906 


0.0015 


0 . O'^go 


0 . 7PO" 7 


O 0 333 


O . 0087 


4.00 


0.3906 


0.0016 


0 . 0780 


O . 3H78 


0.0341 


0 , OOPS 


3.90 


0.3750 


0.0015 


0.0300 


0 . 701 1 


0.09 4 A 


O . 008* 


3.80 


0.3594 


0.001 4 


0.0327 


0 . 3737 


O . O 350 


O . 0088 


3.70 


0.3672 


0.0015 


0.031 1 


0.31 19 


0 . 03A9 


0 . 0091 


3.60 


0.3594 


0.0015 


0 . 0722 


0 . 3972 


0 . 0370 


O . 0097 


3.50 


0.3438 


0.001 5 


0.0346 


O. <4 77 


0.03 7 A 


O . 00*4 


3.40 


0.3477 


0.0015 


0.0340 


0.341 1. 


0.0788 


0.0097 


3.30 


0.3398 


0.0015 


0.0359 


0.35 <7 


0 . O W7 


O . 0099 


3.20 


0.3281 


0.0015 


0.0372 


0 . 3737 


0 . 040 A 


0.0109 


3.10 


0.3281 


0.0015 


A . 0372 


0.3737 


0.0419 


0.0105 


3.00 


0.3125 


0.0015 


0 , 0401 


0.402R 


0.0423 


O.01C7 
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I 



3*00 


0*3125 


0*0015 


0.0401 


0*4028 


0.0428 


0.0107 


2*90 


0*3125 


0*0015 


0 * 040 1 


0*4028 


0 * 0443 


0.011 1 


2*80 


0*3125 


0*0016 


0*0401 


0*4028 


0*0459 


0 . 01 15 


2*70 


0*2969 


0*0015 


0*0433 


0*4349 


0*0470 


0*0118 


2*60 


0*2969 


0*0016 


0*0433 


0*4350 


0*0488 


0.0122 


2*50 


0*2969 


0*0017 


0*0433 


0*4350 


0.0509 


0.0127 


2 * 40 


0*2891 


0*0017 


0*0451 


0*4524 


0.0526 


0 . 0 1 . 31 . 


2*30 


0*2813 


0*0017 


0*0469 


0.4709 


0*0546 


0.0136 


2*20 


0*2813 


0*0017 


0*0469 


0*4708 


0.0570 


0 . 01 4 ? 


2*10 


0*2734 


0*0018 


0.0488 


0*4902 


0.0594 


0.0148 


2*00 


0*2734 


0*0018 


0*0488 


0*4903 


0.0624 


0 * 01.56 


1*90 


0*2656 


0*0019 


0.0509 


0*5109 


0.0653 


0,0163 


M- 

CO 

o 


0*2656 


0 * 0020 ’ ' 


0*0509 


0*5110 


0.0689 


0 . 017 ? 


1*70 


0*2656 


0*0021 


0 * 0509 


0.5111 


0.0730 


0.0182 


1*60 


0*2656 


0*0022 


0*0509 


0*5112 


0.0775 


0.0193 


1*50 


0*2656 


0*0024 


0*0509 


0*5114 


0.0977 


0*0206 


1*40 


0*2734 


0*0026 


0*0489 


0*4911 


0.0890 


0.0222 


1*30 


0*2813 


0*0029 


0*0470 


0*4719 


0*0964 


0.0740 


1*20 


0*2969 


0*0034 


0*0435 


0*4365 


0.1056 


0.0263 


1*10 . 


0*3125 


0*0040 


0*0403 


0*4046 


0 * 1. 1 65 


0 . 0^90 
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17. Computer Program 

In the section the FORTRAN program used is listed. 
Results of cur calculations are given in Ref. 6. 



C TWO 




t DATA FILii 



Ul i 


1 ) = 


15 30. 


'JC l 


2) = 


1 J 7 >.3 


U ( 1 


3) = 


9 JO . 


U( 1 


4> = 


9)1.7 


■J( 1 


5) = 


■+• i. . 4 


Ul 1 


i) = 


6 76.3 


U( 1 


n = 


t 7 . 1 7 


Ul 1 


3) = 


9 40 0 . 


Ul 1 


J ) = 


1 . 


U( 1 


1U) 


= 13. 


J( i 


11 ) 


= 1.3 


Ul 1 


12 ) 


= 9 4 9 0 . 


•J( 2 


1) = 


1 5 JO . 


.Jl 2 


2) = 


1039.4 


'Jl 2 


3 ) - 


3 00. 


U( 2 


4) = 


47 7.9 


Jl 2 


5) = 


5 2*0 


Ul 2 


0) = 


7 3o . 5 


Jl 2 


7) = 


4 7 . 3 5 


Ul 2 


3) = 


9 *+40. 


Ul 2 


->) = 


1 . 


Ul 2 


10) 


= 13. 


Ul 2 


11 ) 


= 1.3 


Ul 2 


12 ) 


=04 *0. 


Ul 3 


1 ) = 


1500. 


Ul 3 


2) = 


1 lOo . 2 


Ul i 


3 ) - 


70 0. 


Jl 3 


4) = 


062.3 


Ul 3 


5) = 


o4 • 3 


ul 3 


6 ) = 


7 33.5 


Ul 3 


7 ) = 


4 3.73 


Ul 3 


3) = 


0400. 


Jl 3 


7 ) - 


1 . 


Ul 3 


10) 


= 13. 


U( -> 


1 1 ) 


= 1.5 


U< 3 


12) 


=9490. 


Ul H 


1) = 


: 1 50u . 


Ul h 


< l ) = 


'1115.6 


Jl A 


3) = 


'SCO. 


U ( J 


4) - 


' »4 3 . 9 


-J l 4 


3 ) = 


7 3.4 


Jl ^ 


6 ) = 


'773.9 


Ul 4 


7 ) - 


'40. 75 


J(*t 


o ) = 


'^440. 



KMX, 
u 1 5) 
GMX, 
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uc 


4 


,91 = 


1 . 


U( 


4 


,10) 


= 13. 


U ( 


*+ 


, 11 ) 


= 1.3 


U( 


l t 


, 12 ) 


= 9 0 . 


J( 


0 


, i J- 


1 ft 0 o . 


«J{ 


? 


, t) = 


1133.3 


IJ( 




,3) = 


3 00. 


0( 


b 


, 4 ) •= 


9 26. 7 


U( 3 


, 5) = 


96.7 


U( 


0 


, c 5 = 


-313.9 


J( 


5 


, 7) = 


ft J . /6 


J( 


5 


, d) = 


9990. 


U( 


5 


, 4) = 


1 . 


J( 


'J 


, 10) 


= la. 


u( 


5 


, 11 ) 


= 1.3 


0( 


0 


, U) 


=9490. 


J( 


0 


, 1) = 


1300. 


U( 


u 


,2) = 


1 94ft . 


U( 


CO 


, 3) = 


9 00. 


J ( 


o 


,4) = 


1300. 


u( 


6 


, 5 ) = 


3. 


U( 


o 


,6) = 


12ft. 3 


U( 6 


, 7) = 


74 J . 


U( 


t> 


, fc) = 


4000 . 


J( 


6 


, 9) = 


2. 


U( 


D 


, 10) 


=2.016 


U( 


6 


, 11 ) 


= 1.4 


U( 


6 


, 12) 


=4000. 


U( 


7 


, l) = 


1 ftOO . 


IJ ( 


7 


, 2) = 


19o. 


U( 


7 


, 3) = 


3 JO . 


J( 


7 


,4) = 


1 7 70 . 


•J( 


7 


, 5) = 


1 l. 


J( 7 


* o ) - 


123.7 


U( 


i 


, 7 ) = 


7 30. 


‘J( 


7 


, d) = 


4 000. 


U( 7 


, 9 ) = 


2 . 


U( 


7 


, 10) 


=2.016 


J( 


7 


, 11 ) 


= 1.4 


U( 


7 


, 12) 


=4000. 


U( 


6 


, 1 ) = 


1500. 


u ( 


3 


, 2) = 


19^3. 


U( 


6 


,3) - 


700. 


U( 


a 


, 4) = 


1 7 30. 


U( 




, 5) = 


13.5 


U( 


a 


, ft ) = 


130.7 


J( 


8 


, 7) ^ 


7 50. 


•J( 


8 


, 8) = 


-,000 . 


U( 


8 


, 9 ) = 


2 . 


U( 


J 


, 10) 


=2.016 


U ( 


8 


, 11 ) 


= 1.4 


Ul 




, 12) 


=4000. 


U( 




, 1) = 


1 50 0 . 


01 


9 


,2) = 


2 030. 


U( 


o 

> 


, 3 ) = 


oOO . 



U( 0,4) =1 70 0. 
Ut ' ■ 

U( 

Jt 

u( 

J( 

J( 



*, 5J =1 7.3 
7 » u ) - t 3 3 • 3 
u , 7 ) - 7 5 0 . 
j =4 G'J3 . 
, 1=2 



J f 



2 f 



10i-2. 015 



-», 12 1 =4QU 3 . 
1 u . 1 * - 1 5 o 0 . 

LU,wl — 2ls*0. 



U( 0, 11 1 = 1 .4 
Jt 
ut 

Jt - . 

U( LO ,3> =5JJ 
Ui lu»'t) = 16o3. 
Ut L u , 5) - !■> .5 
Ut 10,61-137.5 
Jl 10,71=750. 

U( 10,3) — 4J00 • 

U ( 10*91=2. 

Ut 10,131=3.016 
‘Jt 10,11) = L .4 
U( 10 , 1^ ) --t-003. 
7 1 1 , 1 ) = 1 • 
7(1,31=1. 

V t 1 , 3 ) - 0 . 

V ( 1*41=0. 

V ( 1,51=1 . 

7 ( 1,0 = 1 . 

V( 1, 7) = 1 ._ 

7(2, 1 ) =. 33 

7 1 2 , 2 ) = . 3 £ 

V( 2, j) =.25 

V ( 2 * 4 1 = . 3 
V( 2, 5) = . *3 
V( 2 ,6) = . 35 
71 2, 7 ) = 1 . 

V ( 3 , 1 ) = • 7 

V ( 3 , 2 ) = . 7 

7 ( -i , 3 1 = 1 • 

7 ( 3 , A ) = 2 . _ 

7 ( J, 3) =• 93 
7 ( 3 , t ) = . 7 
V( J, 7) =1 . 

7 ( 4 , 1 ) = . 9 
7(4,21=1. 
7(4,31=1. 
7(4,41=0. 
7(4,51=1. 

7t 4, o) = . 5 

7(4, 71 = 1. 

1 ) = . 3 
21=1. 
j 1 = 2 . 

i) *• 1 
£. ) = L • 

W( 2,31=5. 



w ( 1 , 
ftt 1, 
w ( L , 
Is ( 2 , 
«( 
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DATA input 
50 0 k R I T c ( 5 * 5 l J ) 

510 F ij K i A 1 ( ' 0 ' f 0 A t ' E 
i^J) i U » l 
52 C FORMAT (312) 

H; ( IU.ifO. 0) 

I F ( I U • u ) .-2 ) 

IF { !U..i J.-3 ) 

IF ( I U . c G .-4 ) 
KRlTcto, 530) 
fu -MAT ( ' O' f ja, 



0 fER 

V » I rt 



CONTROL CODE’) 



52 5 

53 0 

540 



ST JP 

go r.j 

TO 
TO 



GU 

.jO 



130 

300 

232 



R £ aD ( 5 t 5 4 J ) J ‘ J » 2 U 



ENTER FLUID PROPERTY CHANGES, I i 



530 
5a 0 



57 J 
530 



FORMAT ( I 2, F9 


• X- j 






IF ( J G . c 0 . 0 ) 


GG 


TO 


5 50 


IF ( JU.E3 1 ) 


GO 


TO 


500 


IF ( JU.Eo.-2) 


G J 


TJ 


130 


1 F ( J U . c j • — 3 ) 


J 


To 


300 


IF (JO. EG. -4) 


GO 


TO 


232 


U ( IU ,JJ ) = 20 








Uu TO 525 








■n R I T c ( 5 , 5 f> 0 ) 








FORMAT ( * O' , 3X , * E 


;J ft 


R SY 


R!; AO ( 5 , 5 4u ) 


JV.ZV 




I F ( 0 V . E 0 . 0 ) 


GO 


70 


5 70 


I F ( J V . E 3 . - 1 ) 


Go 


TO 


500 


IF( Ji/.cO.-t) 


GU 


TO 


130 


IF ( JV.2 j . — 3 ) 


GO 


T J 


S 00 


IF (JV.E j • - 4 ) 


oJ 


TO 


232 


V ( IV , J J ) = 2V 








GO TO 550 








R l T E ( 6 , 5 o 0 ) 








F w K'l A T l 1 J 1 t J X » 1 £ 


NTE 


K PCI 


* b AD { 5 ? 3-r J ) 


J*tZ 


N 




I F { J 1% • c G • 0 ) 


GO 


TO 


1 30 


I F ( J . . • c J * — 1 ) 


GO 


TO 


500 


lr (JW«cv*« — 2J 


Go 


TO 


130 


IF (\Ja • li w • “* 3 } 


Gl» 


TO 


300 


I r ( J h . t ; • — ) 


GO 


TJ 


232 



PARAMETER CHANGES, 



T CHANGES, IF ANY.') 



( I w , J .> ) = 2 w 
GO TC 570 
13 0 PS ( 1 ) = J< I'J, 1 ) 

T S ( 1 ) = J ( i U , 2 ) 

P( 1) =U( I J , 3 ) 

I ( 1 ) =U ( 10, 4 ) 

GE LH=iJ ( I J , 5 ) 

JE LH i o = J ( IU,o) 

PHC){ 6)='J{ I J , 7 ) 

g a i.i a i = i ! - algg (tsi i ) / m n /aluoi 3 s i i )/ j m > 

j !/i A 1 = 1 . / ij A -IMA 1 
G 1 1= ( Ga.Y •Ul-l.1/2. 

j12 = u1H-1. 

G 1 a — o A *i ^ A 1 / ( 'j A i'l 1 1 A l — 1 . ) 
l-,= 0l3“l. 

j 15 = 012/ ( GAMMA l- 1. ) 



AN Y . ' ) 



IF ANY.') 
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G l rs = G 1 j— 1 . 

CP 1=0ELM/ t I 3( i )-I(l ) ) 

F 1=CP1/ ; li 

.< V = l . c 3 5 F 7 /.< L 

‘<F H J I 0= i r 2 X ( J{ l J , 7 ) ) 

I r ( 'JPlJl J . _ 0 . i ) J.j TO 400 
» • 1 J ( IJ t 1 *_> ) 
uA * 1A<£ = U ( 1 U , 1 1 ) 

C32 = U( i ,, 12) 

GO TC «»10 
400 Oo2=C61 

G A V A A2 - G A 3 M A L 

>‘»0 “ti V 

410 E T Ah = V ( i V, I ) 
cT 40= V i I V , 2 ) 

A3 A + = V l ! V, 3 ) 

C r = V ( I V , 4 I 
ET«X=V( IV, 5) 

E T -» P = V { 17,6) 

P 3 L P S 7 = V ( IV, n 
OcLY = .-. ( I , l ) 

V "110 = '• ( I. -.,2 I 
T ’1 A X — i> ( I V » J ) 

PJ=lA.o9o 

ro=52C. 

*XlTc ( 6 , L-+0 ) Pi ( 1 ) , TS-Il) , °l 1) ,T( l ) , UELH ,0 ELH 1 6 , KriO ( :>) , Cc 1 , 

+ J ~ L ; J 1 C , v u » G A M 1 A 2 ♦ C 3 2 
I-fO FJM -IaTi * J* , 3X, '(-L JIG PPOPcPT I E 3' // 

f-tA,' I PS 1 = 1 , F 7 . 1 , 1 P S I A • , 1 3 X , ' 2 TS1 = , ,F7«1,' DEG R * / 

+ *tX , 1 :> P i = ' , r 7 . I , * P 3 1 A • , 1 a , ' V T 1 = 1 » F X . L , 1 DEG K * / 

HX, 1 5 { ,u L-Hl ) = * , F 7. if • j T U / L j • , 6 X , * 6 { .lSl-Ho> = * ,F7.2 , 

♦ * DTU/ L 3 * / 

+ 4X » ’ 7 Xrio c - 1 » E 1 2 • 4 , 1 LE>/r T3 • ,3 X , • 6 C6 1= 1 , c 1 2 . 4 , 1 M2 DEG K/CM 

+ 3 * // 

<--*X,‘ S N F L 'J I 0= * * 1 2 , 2 1 X , • L 0 *>, = ', F.12. 4/ 

+ 4X , • 11 GAM 3A2= 1 ,r / . 4, io*, • 12 3 4 2 = ' , E 1 2 . 4 , • M2 DEG K/CMB’ ) 

rt K I T E ( 6 » i 5 U ) ETAt,FTA0.-*3A4,CF,FTAX,ETAP,PSlPS7 
150 F JKMAT ( 1 j ’ * 3 X * 1 i Y i 7E l P AKA ME TE KG' // 

+ 4X , ' 1 “I"Ac = *tr3*2iiXf* 2 uTAi)-*,F5#2,3X 

+ • 3 A3/A t= 1 , F 5 # 3 , 2 X , 1 -t C F= • , F 5 . 2 / -+X , • 5 EIAX=',F5.2,3X, 

+ • 6 FT \P= ' ,ri>.2 » JX, • 7 PGI/P>7-' , F 5 . 2 ) 

wr< I T E I •) , 155) DELI , f MIN, VIA A 
155 F JKMAK • O' , 3X, 'C^.-IPurATIu ,AL FCR-iAT*// 

*4X,' I JcLY* ' ,F 5.2, aX, ' 2 YH 1,-1 = ' , F 5.2 , jX , ' 3 Y -1 AX = ' , F 5 . 2 ) 

150 Y = Y AIN 

P(2)=F{ l) 

P( o)=P( 2 ) 

K 2 = 1 .<=-353 7 / W G 
T( 2) =1 ( 1 > 
r ( 6) = T < 1 ) 

G2 1= ( GA -1 3A2- i • ) / 2 . 

32 2- G 2 l 1 . 

G2a= G 1 M A 2 / ( G A M M A 2 ~ 1 . ) 

G2 •+ — g2 j - 1 . 

325= G22/I GAMMA2- l. ) 
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y] 



«26=G25-1. 

a 1( 1)— ( T Si 1)/T( 1 ) — L. 1/311 
X K L ) = i -J kT ( X M I) ) 

/ ^ C '■? = X -1(1) ( GAi 1/iL/GA 1 YA2-RL/ 12 * T ( I ) /I ( 2 ) ) 

F 1 1 = Af ‘ ( i) : ‘< l.tul i •*<•!( 1)**2)*M - > 1 3 ) 

GMLaxMl i ) / i 1 . 1) 2 ) -30% T ( l.*Gll*X.M( 1) **Z) 

ok i r- ( j , 2 j 3 ) / 1 jA^' 1 mIjY 2 Ck 

23 J F j .< t m T 1 • j » , ^ X i * m V = • ,212.4, ' ( £ F F MOL WT VAPQfl) GAMMA 1 = ' , F7 . 4 , 
+ CR = ' » F 7 m 't ) 
u 1 J To p 0 0 

C MC.MENUJM FUNCTIONS 
aoo Y = Y MAX 
6 Ip nRl Tt( o, 32 J) Y 
62 0 FOk.-UT ( • 3* , 4X, 'Y = * , F7.3) 

X— 

830 DO 850 L 20? - 1 » 2 1 
X Sr* f 1 = c X C 2 S S l X ) 
w%ITtE(o,*i40) X , X :> .M T M 

640 K>.iAT(9X» • <= • rF7.3, • XSM TM= • ,212.4) 

X = X + . 05 

85 0 C G i\ T 1 N j E 

IF ( Y .Lc. YM IN) i3U TO 500 
Y= Y-DE l Y 
GO T G 0 1 U 

C tou ILlrJKI'J -1 AND PERFORMANCE 
2i<: wRIT2(5,2p3) 

23 5 FORMAT i • J • , 3 < , • Y ' , iX , ' X' , 5X, ' A i/Ai • , 3X, • OVER AL L' , 

+ 2 a » 1 RELATI vc • , 3X , ' 13 ' ,oX , • M5 * / 

+-3/X, 'ErF' , 7 X , 'EFF' ) 

Y = Y A X 

20 0 >=XZERU(Y) 

300 P 3 7 = 3 S i 1 ) /P jIP37 
v.c = C?3- JcLT 3 

*P =1^4./ 778. 2 6*1 FS7-P(6) )/RH0(6) 

*EX=«fc/X 

IF ( aLX.LE.iOELh) GO TO 3^0 

WEX-OttH 

aRITE (6,320) 

320 FJk-UT 1 'o' ,*-X, » AVAILABLE ENTHALPY DROP ATTAINED' ) 
j /, J E T - E T A X - « E X - v- P / £ F A P 
2IN*ufcLHlo-rfP/ETAP 
E TAG -'..NET/ jlN 
a T AX = fc T A X *OcLH- ^ P/E T AP 
£ T A R = v. J L T / f, -1 A X 

.•.k[TE ( p , 3 p 0 ) Y, X , U A3»ETAC ,ET \k ,X'-U 3 ) , XM< 5 ) 

35J FORMAT! ' 0' ,Fv.2 , pF10.4) 

Y= Y-OELY 

IFIY.Gl.Y i I N ) vd T G 238 

GO T C. 5 0 0 

END 

C F J vC T I GN XMACH 

fjncT1.,n xnaghi f in, g » 

IMPLICIT KEAL-41 m-.(, u-Z) 

MFLAG-J 
XGUN T = 0 
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JO TO 10 00 



,5,* COOKING L I '.IT EXCEEJEOM 



if- TER 20 I TER AT IONS’ ) 



ul=< J-i.3/2. 
o2-G 1 + L . 
jJ = J2/(J-1. 3 
= 5— L • 
a.-1N= 0 . 

F M A X — G 5 ) 

r- -1 w r v X = F -lx / F M A X 
II- ( F 1C.- iX.uT. 1. ) 

X -1 AC 1 1— 1 . 

«RITC(Of 10 I 0) F -INF -IX 
L 01 0 F ] r '-'ATI • , JX, • F.-iN/r-MAX=’ , F 7 

RE TURK' 

LOGO F=F--N\*l I .+0 L*a.-.N**2 ) **G5-XMN 

IF ( MFLA I.tO. i) . w k I T C ( j i 1 D 2 0 ) XM\,F 
L 03 0 FOR-IAH • J* , 5 X , * -I = ' ,F 7.5, » F=’,E12.4) 

Jb lF = A d S ( F j 

IF { JCLF . IT . . _>E-05) GU TO 1020 
F ^ — FMk •* 2 * X -i ! { L • +3 i *XM\ 2) ^ 36- l . 

X Ah= <M ,-F/F F 
<JUNT=,\ JUNT + l 

IF ( KLUi.T .Ll. ?J) JJ TQ 1000 
aUTEIo, L KO) X *1 iN » F 

10'tO FJrHATt '3' ,<+a,M = ',F7.5, • F = ',E12.4, 

1020 XM AC H= X 'A N 
RE TUinN 
FmD 

; FUNCTION XZcFJ(Y) 

F JNC TU.j XZEROm 
HPLICIT acAL-4 (A-H.O-Z) 

I T E R = 0 
XR = 1 . 

E= cX CE .> $ ( L . 3 
A'3SE = AJS(E) 

I F ( A 0 S £ . 3 E . . 5 E -0 4 ) JO TO 700 

XZckOI. 

RE TURN 

IFIE.GT.O.) GO TO 720 
»RITb( 5 , 71 o ) b 

710 FJkRATI • O’ , 4X , ’ X=l. 0 E= ’ » EL2 . 4. ’ MOMENTUM ObrluIT*) 
<ZERC=l. 
kE TUhN 
720 XL = 0 . 

c - rXtc jS ( 0 . ) 

Ad S E= A3 S ( E ) 

IF(Ad^E.Gc..5b~0‘t) GO TO 73 0 
X2 EP.O= 0 . 

AC TUruN 

IF ( 6 .LT. J. 3 GO TO 7 50 
R I T c ( o , 7 A j ) E 



70 0 



/i J 



740 FUinMAH ‘ J’ , 4X, ’X-0.0 
aZ EkC = J . 

RETURN 

75 o a=( XL<-aR)/2. 

E= ExCES SIX) 

AQSE=Au S { 6 J 



E= * . El 2 . -r» ’ MUMEN TUM INVERSION’) 
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1 1 J 



1 o U 



76 0 



770 



1 [EP= ITER + l 
IF ( ITciN.Lt.li) 
<fi l T E ( O « 7 JO i < i E 



GO 



-t A t 



- i 



X - 



TO 7d0 

F 7. 5, ' 



. 5t-0 4) GO TO /bO 



) GO TO 770 



i ■ j .a , • i A T ( 'o' 

<Z ERC=X 
•>c TURN 
1 F ( AGSi. .GE 
XZtkC-X 

aE turn 

iFit.ir.u. 

XR = X 

GO TO 7 5 0 
XL = X 

GO TO 75 0 
EiM D 

t-J ACTION E < C c 3 S ( X I 

FJNCTI J l E aCESSI X ) 
l ILLICIT *EAL * 6 I A— H 1 0-2. 

COM M f i . J iic),TS(o),XM(o),? 

+ U A 3 , A 3 , CP 1 . C P 2 , C P i , C a I 

^uuLr^fjcLIS^iwtcl^ct :T« 0 » c i AX 
+■ F -1 Y » F i'i A X » !- j'l 3 F M X > F 5 "> F u > G A ‘ i .1 A 1 » 



c = ‘ f E L A 



'* t 



AFTER L8 ITERATIONS') 



) 



| , T{ o) , 7 HOI o},S(u) » V ( 4 » 7 ) » 
itbZ » 3 » C F f.)cLH»0c.LHL6»0£LYt 

’ . E TAG * r v ' 1 » F i 2 » r i 3 » F M4» F M 5 » F MX , 
i A 2 »(>’>••! MA 3 i j1 L , ji.2»313»Gl4»Gl K - 



ET 
■j A 



L6 i G 2 L f 'j 2 2 » 23 » o 2 » » j 2 5 » G2 0 » ■ j i * » .*3 2 t G-i 3 > j3 *+ t ->3 5 > l j 3 o » GO 1 » G i 2 t G M X » 
+■ a -I X t X ' '1 Y t .K io » a K 1 » ' r i_ ^ I 0 * P V T 6 » P '.) » P S 1 P S 7 » ? I I » I » a 2 * R 3 » 

G.i I X » T 0 » i nE X , W i ; i *wc f » WG , WV , V , 7. -II N , V MAX , T L , XM L , Jl , FMZ , ?l » 

+ JI 1.) , 12) 2, 3) 

\.i (2 ( X '.( L ) / Y) '-.SORT! { G A IMA 1/ GAMMA 2 )*<Rl/*2)*(T(L)/T(2))) ' 

IF (XM{ 23 . L7. 1. ) GJ TO 2 3:) 
a -I ( 2 ) = L . 

HZI-Tll )/Y**2*^4.-HAL/GAM-IA2*Rl/R2 
23 i TS2T 2= L . +32 L-X.MI 2) **2 
T o I 2 ) = T(2)-TiiT2 
Pb(2)= P (2 ) - TS2T2**G23 

Gi2=X,"I 2 ) / I I . +GA,-!MA2 -a.H 2) --2) *53RT( l.+G21*XM( 2)**2) 

240 «3=A^ t + U.-.<l*iva 
CP2=G23*R2 

CP3=X*0P IK I.-XJ *C?2 
GAXMA3-CP3/ <CP3-a3) 

G3 1= I GA.-V-1A3- I . ) / 2 . 

G32=G3I+I. 

G33 = GA. i VA3/I GAM4A3-1 . ) 

G3 4=G33- I . 
y 3 5— G3 2 / I GA A !A3-l.) 

G 3 u= 3 3 ‘3— L . 

F i A X — 'jJ2^' x {“u35) 

TSI3)=t<*CPL*TS{l)+(L.-X)*CP2*TSV2) )/CP3 
FM 2 = XS( 2 ) •? ( I,+j21 *aM ( 2) -s=*2 > **{ -G2 5 ) 

01 = P S(2)/Pi( l ) -* F v 2 / r , v . L * $ 0 R T I RL/R2^T$( i) / TSI2) 

AI A3 = X$0 1/ I L .-X+X*OL ) 

SI I J=C PLv“;lO ^ (IS i 1) / rj)-Rl*A LOG! PS ( l)/=»0) 

SI ?.) =CP2 -AlO'-I TS (2) /T J)-a2*AL0G(PS< 2 )/P0) 

P3X=(lP )6SLJ j ( To I 3) /TO)-a-S I 1) -( L .-X ) *3< 2) )/ 73 
PSX=PG*F. <*M PSX) 

0 2 = P 6 I 2 ) /PiX«FM2*SiKT( ,\3 /< 2* T S I ) / T S I 2 ) *GA. A2 /G VA MA 3 ) 

FMX=D2/( L.-X + X-JU 



'GA Vi A 2 / G A M M A 1 ) 
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IF (F VX. L£.f 1 AX J GJ TO 2 *L 
X M A = l . 

t- -1 X F > v > = r v A / F M A X 
FMX = F,MAX 

•’Is 1 T t l <> » 2 4 3 ) Y * x f r X XF MX 

2*p FU.V.:«T ( • j« ,4X, * ,F r.5, • X = * ,f 7. 5 , 1 F MX / F MA X- ’ , £12.4) 

GJ TC 24* 

241 X-1X=XM*:H< F M X » G A , >U3 ) 

24 4 GMYxSwiU (,M4TS( 3 ) / G A 1 -1 A 3 ) / ( X /G 11 *S JR T (ML *T S ( I) /GAM W1)K 1.-XI/0M2* 

+ j'jis T ( K2 * T _' ( 2 ) /b«A i v i -1 A 2 ) ) 

GM AX = 1 . / S Jin T ( 2 . * I G*M -i A 3 <• l . ) ) 
i F ( GMY . L c . s»MAX ) TO 243 

Jl Y o v X=GMY/G max 
V, R i T c l f j » ! Y » X , ji'i YgMX 

242 FORMAT ( ‘ )’ ,4X, ’ Y=* ,F 7.5, ' X = ',F7.5,' GMY /3M4 X= • , E 1 2 . 4/ 

+ PX, 1c. iTUP LIMIT EXCEEDED. • ) 

X M Y = 1 ♦ 

F >•', Y = F M A X 
GJ TC 243 

24 5 X -1Y= { ( 1 . -2 . <=GA.XMA3*GMY**2 ) -SQ2T{ 1 .-2 .*( GAMMA j + I . ! *GMY**2 ) ) / 

+ ( l • — o A 5 M \ 3 f ( l . — 2 • * j n 1 1 A 3 v G H Y * * 2 ) ) 

XM'/= 5 }k T ( v-i / ) 

FMY-a t.< t ( 1 . * j 3 i*AMY* -2 )-*G35 
IF ! FMY . LE. r.-UX) vjJ TO 243 
X M Y = l . 

r.-n r v X = F Y Y / F M A X 
r.M Y — F iM -* < 

A K I T E (5,2461 Y , X , F M Y F IX 

2*6 F Jis 1AT { 1 J ' , 4X , 1 Y- ' , F 7. 5 , ’ X*' ,F7.5»* F MY /F MX = • , S 1 2 . 4) 

24 j FSv = ps A-FM x/rMY 

?G( i) = PS-<*( ( PSY/PSO**U./ETA£)) 

FS3=FSa/k${ 5) -F MX 
IFIFM3.LE.F MaX) Gj TO 244 
X -U 3 ) = I . 

Mirt| J( 9jO) X -1 ( L ) , X *( 2) , < IX ,XMY, XM( 3) , F Ml ,F M2 ,FMX , F MY , FM3 , 

+ F 1 A < iu'L , GM2 , GMY , GMAX , PS (1 ) , PS (2 ) , 2 SX , P SY , P$ ( 3 ) , T J ( I > , 

<-TS(2i ,TS( 3) 

901) FOR M A T { • J* , 6X , • M I , 12 , MX, *'Y , M3- ' /5E 12 . 4/ 

** 0* , ' F 11 ,F M2 ,F-1X ,F MY , FMi , r .M AX = • / 5 E 1 2 . 4 , £ 1 1 .*/ 

+ * (J * , * bt'i i , i 'i 2 , M i t 3 I X = * / 4 c 1 2 • 4 / 

+ 1 J 1 , ,D 3l,Mi2, u SX,FjY »PSi =, /5£l2.4/ 

«•* J* , • T jl , I 22 , TSJ = ' /3EL2.*) 

F 13 F 'X = FM 1/FMAX 
F M 3 = r "imX 

-.KITE ( 5 , 2 3 L ) Y,X,FM3FMX 

25 i t-CisMAT { ' J' , 4X , ' Y= ' , F 7. 5, ' X=',F7.5,' F K3/F MA x= ' , cl 2 . 4) 

GO TC 2-, 7 

24 9 Xi-1 ( 3 ) - a 1 A C n ( F M 3 , G AM M A 3 ) 

2-, 7 T S 5 T 3-= i . + G 3 1 ■* X -1 ( J J * ~ 2 
T(3)-Ti(3»/TSJT3 
PS JP i=TS jTp#‘?G3J 
P ( 3 ) = R 5 i J ) / P S 5 P 3 • 

PS ( 4 ) — P i ( 3 ) *iX? ( -{ i ,-fcTAJ) *GAMMA3/2 . *XM( 31 **2) 

F M4 = P S l } i / PS { ,) ^ A J A 4 « p M 3 _ 

[ r ( F Mh .uc.F’IAX) GJ Tu 252 
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n 



X 3 ( 4 ) = L . 

F.HF v X-Fi34/F. 4 AX 
F,3h.= FpAX 

*4 IT F (o , 253 ) Y,X,F 14FMX 

25 3 F ‘J .< *-• A T ( ' J ' , »X. , Y- , »F7»p»* X = ',F7.5,« FK4/F 3A X= • , £ L 2 . 4) 
oj TU 254 

252 X IUIsa'UwHIFM'm j AM M A 3 ) 

25+ T 3 ( 4 ) = T 3 { 3 ) 

T 3 4T 4= L . +G 3 1 v Xi3 ( 4 J ♦* 2 
T i -+ ) = T ;>( 4 ) / T 34 T4 
P i -t P i= T :> 4 T + * * V7 3 3 
P{ 4 )=Oi( 4) /P54P4 
Co i=X*CG l H l.-X)*C32 

x TA = G. JGJVta75*( pJ/L4.69t>) *< Cc33/9490 .*520. /TO) **2 
DELTS / i- = >5£.TM?Gjl/bA i, v i-.3*PS( 4) /P0*( T0/TS14) ) **2/ ( L . +G3 1 * X.‘M 4 ) ** 2 ) 
•*■ G3 4 

OEL T 5= T o< •+) *JELT54 
T5 l 5 ) = T 3 { 4 ) -GEL T 3 
P 3 { 3 ) = P 3 { «•)•*( T o ( b )/l 3(4) ) * *G 3 3 
F 3 5= TS ( j ) / 7 3 1 4 ) *P S( 4 ) / PS ( 5 ) *FM4 
IF (FM5.L=.rMAX) GO TO 250 
XM ( 5 ) = 1 . 

F>3 5 Fi v X= F -1 5 / F,3AX 
F.3 3 - F .SAX 

WAITE (o,255) Y , X , F X5FMX 

255 FUA-1AT l » 0’ ,-+<,« Y^ F7.5, * X=»,F7.5,' F -15/F MA X= * , El 2.4) 

GO TO 2 o 0 

250 X.1 ( 5 ) = X 3aCH ( F45 ,GAM,YA3) 

2b J TSpT 5=; . + j o 1 * a >• ( 5 ) **2 
T ( 5 I - 1 3 ( 5 ) / T 35 T5 
P 3 5P 5= T 3 5 T 5* ^3 3 3 
P( 5 ) = P 3 l 5)/P3^P5 

F 3 Z = ( 1 .-X) *A5A4*F‘32 / ( l.-X*X*Di ) 
aMZ- xF.mCH ( F 3 Z ,GA.V3A2 ) 

TSZT/ = 1. *-G2 1 *XHZ **2 

TZ=T3(2)/TSZTZ 

PZ = ?3(2 )/(l 3 Z T Z * *G23 ) 

fcxCt SS=1 .-PZ/P( 3 ) ■<•( L .-X-CF/2 . ) *GA 3f-‘A3*XM (5 ) **2-( 1 . -X ) *GAMMA3* 

+ X.3 < 5 ) ^X.3Z*30kT (GAM.3A2/3AM 3Ao*R2/X‘3*TZ/Tl 5) ) 

RETURN 

End 
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18. Nomenclature 



A 


area 


a 


speed of sound 


C B 


characteristic breakdown constant (Eqn 10.3) 


s - 


specific heat at constant pressure 


c f 


friction loss coefficient - condenser 


E b - 


breakdown electric field strength 


E (x) - 


Eqn. 11.5 


F(M n )- 


Eqn. 5.16 


f (M) - 


Eqn. 5.3 


^max 


Eqn. 5.15 


g(M) - 


Eqn . 7.5 


g 


Eqn. 7.7 


g max 


Eqn. 7.9 




Eqn. 7.10 


h 


enthalpy 


I 


irreversibility 


M 


Mach number 


m 


mass flow rate 


?e ’ 


gross electrical power 


P 


pressure 


g in - 


heat input 


R 


gas constant 


R 


universal gas constant 


s 


entropy 


s m 


entropy of mixing 
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T 



temperature 



V 

W 



* 



w 



e 



* 



w 

w 

X 



p 

net 



y 

z 



velocity- 

molecular weight 

gross work out per unit mass 

ideal pump work per unit mass 

useful electrical output per unit mass 

V*3 

v l /v 3 

quality of vapor 



6 - Eqn. 10.7, dimensionless breakdown constant 

y - ratio of specific heats 

A( ) - change of 

e - permittivity of the medium 

n c - overall cycle efficiency 

n„ - Eqn. 8.13, ejector effectiveness 

ht 

n D - Eqn. 9.6, diffuser effectiveness 

rip - pump efficiency 

n x - excitation efficiency 

p - gas density 

ip - steady flow availability 

a) - mass ratio of condensible primary to dry secondary 
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Subscripts 



0 - ambient conditions 

1 - primary, entrance to ejector (exit from boiler/superheater) 

2 - secondary, entrance to ejector (exit from separator) 

3 - mixture, exit from ejector entrance to diffuser 

4 - mixture, exit from diffuser entrance to generator 

5 - mixture, exit from generator entrance to condenser 

6 - primary, exit from condenser entrance to pump 

7 - primary, exit from pump entrance to boiler 

f - saturated liquid 

G - gas or secondary 
g - saturated vapor 
n - next value 
s - stagnation 
v - vapor 

x - exit station of hypothetical minimum loss device 
y - exit station of ideal ejector 
z - exit station of condenser 
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